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We study a scalar-tensor version of Lovelock theory with a non trivial higher order galileon term 
involving the coupling of the Lovelock two tensor with derivatives of the scalar galileon field. For a 
static and spherically symmetric spacetime we extend the Boulware-Deser solution to the presence 
of a Galileon field. The hairy solution has a regular scalar field on the black hole event horizon 
and presents certain self tuning properties for the bulk cosmological constant and the Gauss-Bonnet 
coupling. The combined time and radial dependence of the galileon field permits its horizon regu¬ 
larity. Furthermore in order to investigate the effects of linear time dependence we find spherically 
symmetric solutions in 4 and 5 spacetime dimensions. They are shown to have singular horizons. 
Afar from the Schwarzschild radius and for weak higher dimensional couplings the solutions are 
perturbatively close to GR representing exteriors of GR like star solutions for scalar tensor theories. 

PACS numbers: 


I. INTRODUCTION 


Since a number of years there has been a renewed interest and intense activity in the subject of modified gravity. 
This is true not only at high energy scales but also at infrared scales and very large, cosmological distances. This 
quite recent intense activity was initiated with the advent of braneworlds where our universe is pictured as a brane 
in a higher dimensional bulk spacetime. Ultraviolet modifications to 4 dimensional gravity where elegantly predicted 
in the Randall-Sundrum type of models [1] with a brane localized 4 dimensional graviton whereas large distance 
modifications were implemented by the GRS and the DGP models where the 4 dimensional graviton was only quasi- 
localised becoming of higher dimensional nature at very large distances [2], Although General Relativity with a 
cosmological constant is the most general metric theory 1 in 4 dimensions this is not so in higher dimensional theories 
and implications of this to braneworlds were also investigated (see for example [5]). Implementing 4 dimensional 
effective theories inspired from braneworld models gave a new perspective to 4 dimensional scalar tensor theories 
and also more recently massive gravity and bigravity models [6]. In this context while investigating the cosmological 
constant problem it was noted [12] that the most general scalar tensor theory with second order field equations 2 was 
known and written down explicitely by Horndeski [7] back in 1974. In fact, as we will argue now, Horndeski theory 
is closely related to Lovelock theory [8] and this is the direction we will take in this article. 

Lovelock theory is the most general metric theory in arbitrary dimensions with second order field equations [3]. 
Reducing the theory down to four dimensional spacetimes it gives the most important Horndeski [7], terms especially 
all those involving scalar-tensor interactions [11]. Interestingly all four scalar tensor interaction tensor terms are the 
unique terms of Horndeski theory which self tune in flat spacetime an arbitrary vacuum cosmological constant [12]. 
Furthermore, consistent Kaluza-Klein truncations of Lovelock theory give Galileon black holes [13] from their higher 
dimensional images. The first higher order term of Lovelock theory (after the cosmological constant and the Einstein 
Hilbert term) is the Gauss-Bonnet term, 

Cgb = V—g ( R 2 - 4 R a0 R af} + R aMS R a ^ s ) (1.1) 

which is a topological invariant in 4 dimensions. It is also the unique higher order term giving second order field 
equations in 5 or 6 dimensions. Indeed varying (1.1) with respect to the metric gives us the Lovelock 2-tensor, 


H, lv = - 


-9,,, (i? 2 - 4 R a0 R af} + R a0lS R aplS ) ~ 2 RR pv + iR w> RJ‘ + R Kpp „R Kp - 2R l , n x P RJ ;Xp 


( 1 . 2 ) 


* Electronic address: christos. charmous is@th.u-psud.fr 
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1 For the precise theorem due to Lovelock [3] and a extensive discussion see for example [4] 

2 Again for the precise formulation of the theorem one can see for example [7], [49]. See also interesting generalisations of Horndeski 
theory in [9],[10],[59],[60]. 
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which is divergence free and identically zero in 4 dimensions giving us in this case the first Lovelock identity (for 
extensions see for example [14-16]. 

In four dimensions the action involving the Gauss Bonnet invariant when coupled with a scalar field is no longer 
a topological invariant. It has a non trivial variation with respect to the scalar field. It is one of the Fab 4 terms, 
only self tuning with a little help from his friends [12]. It is also a very particular term in the Horndeski action as it 
gives an effective mass term (weighed by spacetime curvature) in the scalar field equation and can bifurcate Galileon 
no hair arguments [17], in the case when the scalar coupling is simply linear for then it still retains translational 
symmetry [18]. Early work adopting perturbative analysis with dilatonic couplings 3 were performed in [19-23], while 
numerical investigations were performed in [24]. Recent investigations on the subject were carried out in [25, 26]. A 
different way to bifurcate no hair theorems is to involve scalar tensor interactions involving translational invariant 
Galileons such as the John term of Fab 4 which reads, G M v S7^(j)S7 u <j) where G „„ is the 4 dimensional Einstein tensor. 
There when one considers additionally a linear time dependence of the scalar field [27] (see also [28-34], and [35-41]) 
it was shown that the time dependence yielded analytic GR like black holes with additionally well defined scalars on 
the black hole horizon. For a recent review on black holes and scalar fields see [42] and [43]. 

In higher dimensions the Gauss Bonnet term as we mentioned earlier is no longer a topological invariant, it is part 
of the metric theory, allowing for static solutions [44-47],[4, 48]. As a result, one can construct other than the 4 scalar 
tensor interactions existing in four dimensions. They have been discussed in the literature but they have not been 
investigated on the type of solutions or effects they may have [49] (see also [50-54] ). Indeed the ’’simplest” of terms 
can be obtained in a similar way as the John term in Fab 4. It boils down to couple the Lovelock tensor H with 
the scalar tensor V M </>V"< j>. The divergence freedom of guarantees that this is a Galileon term (in other words 
yields second order equations of motion) and it reads, 

Co.HGB = y/=g (1.3) 

This term has just like Fab 4 terms [12], self-tuning properties although in higher dimensions the application of hiding 
a bulk cosmological constant is not immediate. Nevertheless it can hide a bulk cosmological constant and we will see 
this explicitly happening here. 

Our aim in this paper is two fold: first we wish to investigate Lovelock theory including higher order Galileons 
as (1.3) and finding when possible scalar-tensor black hole solutions. We will accomplish this using the construction 
techniques developed recently [27], [30],[33]. Secondly we want to study in more depth the effect of time dependent 
linear scalar fields. Towards this aim we will consider the extreme situation where the scalar field has only time 
dependence. Here, the scalar or axion is similar to the solutions found recently for three forms in AdS spacetimes in a 
different physical context more akin to string theory. Can we have regular black holes with a purely time dependent 
axion? Will there be a black hole horizon as for the axionic solutions of [56], [57] where the axions are switched in 
the space-like directions of the horizon? The aforementioned axion solutions have been discussed in the context of 
holographic systems as describing momentum dissipation in the boundary of AdS [58]. In fact the free axionic scalars 
break translational invariance in the horizon coordinates and as a result momentum in the dual three dimensional 
field theory is no longer conserved. The time dependent axion solutions would then in complete analogy describe 
energy dissipation for the dual field theory. We will demonstrate in fact that one has to consider also non trivial 
radial dependence [27] in order to have black hole solutions. We will find that although the solutions have actually 
a well defined GR limit the scalar tensor solutions are singular, unlike the case of the axionic black holes found 
previously [56],[57]. We will also show that the solutions are valid GR like star solutions as long as we are far from 
the Schwarschild radius and for weak higher dimensional coupling. 

Keeping these objectives in mind we will consider the following theory, 

S = j dx d ^j[(R - 2A + a (i? 2 - 4 R a pR af) + R a p jS R a ^ 5 ) - , (1.4) 

Notice that the scalar field has translational invariance. So here, we are not in the context of string theory or 
Kaluza Klein actions where </> dependence is explicit and the translational invariance is broken. This symmetry will 
be important for the partial integrability of the system. We can also see that the dimension of the scalar field is 
[</>] = For the rest of the coupling constants we have that [/3] = —2, [J] = —4 and a is the usual Gauss-Bonnet 
coupling with [a] = d — 4. This is a rather involved gravitational action but we will see that in order to construct a 
scalar tensor version of a Boulware-Deser solution all above couplings are necessary. 


3 This analysis includes the case of linear coupling 
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Varying the action (1.4) with respect to the metric we have the following equations of motion 
CG t] TjSpviPp) ^ 4~ fZ/u' A 4~ 

+/? (hopfc^ + P mv pS7 a ^<j> + 


+s( - \a^H a ^j>V 0 <j> + H { , p \7 v) (j)\7 p cj) - 

-^ClS^^R^Ra, p ^V p <fiV x cj> - ^V^V^Va*) = 0, (1.5) 

where P a p p v is the double dual of the Riemann tensor, P a p P i, = -^\^ a ppaR pai& Cpu-yS, which can equivalently be 
written as P a/3 Rj^.. Expanding the previous expression we have that 

Pappv = Rafipv “4 RppQaii PpuQap Raptlfiv 4 “ RavQpp. 4 “ ~^PdapQpif ~^PQPixQolv ( 1 - 6 ) 

The variation of the action with respect to (j) can be rewritten in the form of a conserved current and this is a 
consequence of the shift symmetry of the action, 

V M J p = 0, J p = (rjg pu - 0G p " - 5H pi/ ) d v c\>. (1.7) 

Note in the above current the presence of the Einstein-Gauss-Bonnet plus cosmological constant metric equations of 
motion. This is a critical condition in order to find analytic solutions of the scalar tensor system. 

In the next section we will consider purely time dependent scalar fields in 4 and 5 spacetime dimensions. In 4 
dimensions the above action boils down to a rather simpler one since the Gauss Bonnet and the Lovelock tensor do 
not yield any part in the field equations, i.e. a = 5 = 0. After finding several solutions we will argue that these can 
never be black hole solutions unlike the case of magnetic axions [57]. We will then go ahead to 5 dimensional scalar 
tensor Lovelock theory and there find explicit black holes with a non trivial and time dependent scalar generalizing 
the celebrated Boulware Deser solution [44]. We will summarize our results in the last section. 

II. A TIME DEPENDENT AXION AND THE GENERIC IRREGULARITY OF THE SOLUTIONS 


We want to examine the existence of spherically symmetric static metric solutions where there is only a linear 
time-dependence in the scalar field. This is a case akin to the electric axionic black holes found recently in [56],[57] 
or to cosmological galileons where linear dependance is customary (see for example [55]. The Galileon in this case is 
an axion switched on in the Killing time direction. It is easy to see that in 4 dimensions its dual three form is in fact 
magnetic. Starting from 4 dimensions the action (1.4) reduces to, 


S = J dx 4 ^[CR^2A-r](d(l)) 2 + pG piy \/p(/)\/ v 


In [27] a linear time-dependent and radial scalar field was considered and regular black hole solutions were found. 
Here we consider, 

ds 2 = —h(r)dt 2 + ——dr 2 + r 2 dfl \ (2.2) 

f(r) 

4>(t,r) = mt. (2-3) 

Doing so the scalar field equation is automatically satisfied. Additionally there is no off-diagonal term in the metric 
field equations. Hence we are only left with the (rr) and (tt) equations for our two remaining variables to be 
determined, / and h. The (rr) equation reads, 

rfh' _ 2 2hA — m 2 rj 2(h + pm 2 

h 2 (h — Pm 2 2 (h — pm 2 


whereas the (tt) equation reads, 


V(f - !)]' = ~r 


2 2hA + m 2 g 
2 (h — pm 2 


(2.5) 
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In the following we will examine various cases for the values of the couplings, always supposing that the coupling 
constant p is non-vanishing. The system of equations is not integrable but we can make some qualitative remarks. 
We remind the reader that the ( coefficient denotes the ”GR” - part of the field equations. Then as we can see the 
effect of the Galileon term is to ’’shift” the zero’s of h. This as we will see leads to singular spacetimes as possible 
zero’s of / are not zero‘s of h. 


A. Four dimensional solutions 


Let us start by 


r) = A = 0 and ( ^ 0. 

which is integrable. In this case from the (tt) equation we get, 

(2C h-m 2 p) (—1 + / + rf') = 0. 

2 o 

If h(r) = ^2- then we need to have /(r) = 1. When the second bracket is zero we have 

(-1 +f + rf) = 0 /(r) = 1 — —. 

r 

On the other hand the (rr) equation can be written as, 

[ u — (3m 2 _ l u \ 

J Cfl (2(h + p m 2 )h J d; f 

which can be easily integrated and solved for h. We thus have, 

(2 C,h + to 2 /3) 2 = cfh , 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


where c is an integration constant that fixes the asymptotics of the solution. Clearly as / = 1 for large r we need to 
fix c = (2£ + m 2 / 3) 2 so that h = 1 for large r. The solution for h reads (only the upper branch is valid), 


cf — 4m 2 /3£ + vW cf — 8 m 2 /3( 

8C 2 


( 2 . 11 ) 


for the above value of c. Here we note that taking pm 2 —>• 0 the solution smoothly flows to the GR Schwarzscliild 
solution. Hence we are in the same branch as for GR. This is unlike the stealth black hole found in [27] where the 
branch is distinct (although the metric is identical to GR!). However although the solution here is connected to the 

2 a 

GR solution it is not a black hole for when we hit the event horizon at / = 0 we have h = The solution is 

therefore singular since It will always develop a branch singularity at the zero of the square root when / < 0. But as 
long as we are far from the Schwarzschild radius, as is the case for a spherically symmetric star and << 1 the 
solution will look like an exterior GR star solution. 

The generic case where all coupling constants are not zero is not integrable. However note from the (tt) equation 
that setting ^ + A = 0 yields a common root on the RHS giving immediately, 


,/ x -i M , 1 2 

f(r) = 1 - - + —r . 
r 3p 


Now we can proceed as in the previous case since the (rr) equation reads, 


and can be written again as 


I" 2ph — Pm 2 
J (2 C,h + Pm 2 )h 



/ + ^> 



2 (,h — Pm 2 
(2C,h + pm 2 )h 



( 2 . 12 ) 


(2.13) 


(2.14) 
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The solution for h(r) comes from solving the algebraic equation 

(2 (h + m 2 /3) 2 


= c/, 


where c is an integration constant, or 


h = 


cf — 4m 2 /3£ ± vW c f ~ 8m 2 /3C 
8C 2 ' 


(2.15) 


(2.16) 


Again we conclude as before that the solutions are not black holes, although again they have a smooth GR limit. 

Since we are not able to analytically solve the field equations for the full set of coupling constants it is instructive 
to adopt a perturbative approach in order to find the leading order behavior of the solutions. This is going to be 

2 a 

performed by expanding the metric functions in terms of the dimensionless parameter e = supposing that e << 1 
and solving order by order the field equations. Additionally we will suppose for simplicity that 77 = 0. We will thus 
have 


f(r) = fo(r) + e/i(r) + 0(e 2 ) and h(r) = h 0 (r) + ehi(r) + 0(e 2 ) where e = 


m 2 f3 


«1 


(2.17) 


From the zeroth order term we find that fo(r) = ho(r) = 1 — (r — resulting in an (A)dS spacetime, as expected. 
Going to first order we find that 


= - 

hi(r 


Cfo{r) 

f fi(r) 


fo(r)J \fo(r) 


2 (/o(r)-l) 


rfo{r) 2 


Integrating the right hand side of (2.18) we find 

/' 3r 3 A 


= 3 r + 9C 


(.Xi — to) ln(r — Xi) 


J -3rC + r 3 A + 3 (fi 
where Xi are the roots of the polynomial 

Ax 3 — 3(x + 3 C/U = 0. 


3Aa:| — 3<( 


(2.18) 

(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 


This is sufficient to see that our approximation breaks down when we approach any zero of fo since /1 becomes clearly 
dominant over fo close to the horizon. Therefore we come to the same conclusion as we did for our specific analytic 
examples. The solutions have singular horizons. 


B. Five dimensional case 


For completeness we will examine the case of purely time-dependent scalar field for the full action (1.4). We start 
by choosing a five dimensional static metric of the form 


ds 2 


—h(r)dt 2 


1 

f(r) 


dr 2 +r 2 dnf 3U , 


( 2 . 22 ) 


where dfl 2 represents a three-dimensional maximally symmetric space, where k = — 1 , 0,1 is the spatial curvature 
of the three-dimensional hypersurfaces, while for the scalar field we have 


<f>(t, r) = mt. 


(2.23) 


The (tt) equation reads, 

—6(to 2 (5 — 2ah)[(2n — /)/]' + 2r(3«;(m 2 /3 — 2(h) + r 2 (m 2 ij + 2 Ah)) — 3(m 2 /3 — 2 (h)[r 2 f]' = 0 


(2.24) 
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and the (rr) equation reads, 

— 12 (k — f)f(m 2 S — 2 ah)h' + 2 rh( — 3(k — f)(m 2 (3 + 2(h) + r 2 (2Ah — to 2 ??)) — 3 r 2 f(m 2 (3 — 2(h)h' = 0. (2.25) 

It is easy to see that setting a = S = 0 we end up in the higher dimensional version of (2.4) and (2.5). We now have 
the following constants A, (, a, rj, (3, 6 and k for the spatial curvature. Setting all coupling constants to zero apart from 
6 leads to trivial solutions for f(r) and h(r), (mind also that if f(r) = n, then h(r) is undetermined). Let’s examine 
the following case. 

We set: 


(r] + (3 A = 0 and ar] + 6A = 0. 

The ( tt ) equation gives, 

(to 2 ?7 + 2A/i)(2?’ 3 ?7 + 6 r(3n — 6 r(3f — 3 (r 2 /3 + 4 8k — A8f)f) = 0 
while the (rr) equation becomes 

2r(r 2 r] + 3(3k — 3(3f) (?n 2 ?? — 2Ah) h + 3 f(r 2 (3 + 4 6k — 4<5/)(?n 2 ?? + 2A h)h! = 0. 
The solution for f(r) is, 


(2.26) 


(2.27) 


(2.28) 


f(r)=K+ — ± 


r 2 (3 r 4 (3(3 2 — 4 6rj) 


48 6 2 


fo- 


whereas as before h is solved by the quadratic, 


(2.29) 


(h-^) 2 = cfh (2.30) 

and as before we can conclude that there can never be a horizon covering any of the branch or central singularity 
at r = 0. It is absolutely essential that the scalar field has to have a radial profile. This not only regularizes the 
geometry but also the horizon behavior of the scalar. We move onto this now for the case of Lovelock galileons in 
D = 5 dimensions. 


III. SCALAR TENSOR LOVELOCK BLACK HOLE 


After examining the axion time dependence on static scacetinres we will now proceed to the case where the scalar 
field has additionally radial dependence. As we will see the combined time and space dependence of the scalar will 
regularize it on the horizon. We will look for black holes to the full Lovelock Galileon action (1.4) in 5 dimensional 
spacetime. 

We choose the same static metric Anzatz as before 

ds 2 = - h(r)dt 2 + dr 2 + r 2 dn%, , (3.1) 

f(r) v ,K 

We suppose that the scalar field reads, 4> = <j)(t,r ) ie., it does not have the same symmetries as the spacetime metric. 
For the above ansatz the only non-vanishing current components of are the following: 


Jt = 2^( Mr 2 V + 3(3K-3f3f)-3(r 2 f3 + 46K-45f)f ) (3.2) 

J r = ^-{2r(r 2 r 1 + 3(3K-3(3f)h-3f(r 2 (3 + A5K-A5f)h'), (3.3) 

2 r A h 

where (') and (') represent differentiation with respect to t and r respectively. We see that, due to translation 
invariance, the scalar field equation can be written as 


Vp 


_ dJ l dJ r 1 / 6 

dt dr 2 \r 


r 

f 



= 0 . 


(3.4) 
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For a generic scalar we see that the scalar field equation can not be written as a first integral and we also end up 
with a PDE system, which even in the simpler case of [27], is not integrable [34], Only when the scalar has a linear 
time dependence, we have a system of ODE s with no explicit time derivatives (but in the presence of the velocity 
charge q ). In fact one can show that this is the only case where we have consistency and integrability of the system 
of ODE’s [33]. The energy momentum tensor associated to the scalar field then does admit a time like Killing vector 
as it should. We therefor have, 


<t>(t,r) = qt + ip(r). (3.5) 

Due to the linear time-dependence of the scalar field, there is still an off-diagonal part (tr) in the metric field equations 
which is not however time dependent. As was proven in [33] for (3.5) setting this off-diagonal part to zero is equivalent 
to the metric satisfying the following equation 

r? g rr - (3G rr - SH rr = 0. (3.6) 

The above, (3.6), immediately satisfies the scalar field equation (1.7) as J r = 0, [27],[30],[33]. Therefore (3.5) solves 2 
of the field equations at the same time rendering the whole system mathematically consistent. This means that (3.6) 
can be solved for f(r) finding that, 


./O’) = 


6 rfPh + 3 (r 2 /3 + 45k) h' ± yj— 96r<5 (3/3 k + r 2 r]) h h' + 9 (2 r/3h + (r 2 /3 + 4 5k) h') 2 


245h' 


(3.7) 


We can now use this expression in order to find a relation between the scalar field function ip and the function h. 
Indeed if we substitute (3.5) and (3.7) to the (rr) component of (1.5) we can solve for ip' . The expression is highly 
non trivial involving a higher order algebraic equation for ip'. But it can still be done and in principle we can replace 
ip 1 in the (ft) equation obtaining a non linear ODE for h which can be eventually solved numerically in general. 

Here to go further analytically we choose to set f{r) = h(r) in the metric and search for solutions in the manner of 
the self-tuning de Sitter Sclrwarszchild solutions found in [27]. From (3.6), it is easy to see that the solution reads, 


f(r) 


r B 1 

'‘< r > = ' t + TT ± I 


r 4 (3/3 2 — 45 rj) 


3 5 2 


16ci. 


(3.8) 


Now from the (rr) component of (1.5) we can solve for ip 1 , substitute it back to the equations (alongside with ip") 
and see that the field equations are satisfied, if and only if the coupling constants satisfy the following relations: 


5(q 2 (3 + (k) — a p) k = 0 
t)((k — q 2 p) + A(3k = 0 


as well as, 


ip'(r) = ± 


g\/(«- /)/« 
/ 


(3.9) 

(3.10) 


(3.11) 


where we have used (3.6) and (3.9),(3.10) in order to solve for a and A. It is obvious that eliminating k from (3.9) 
and (3.10) we get ij(a/3 — 2 5£) — f35A = 0. 

This solution is very similar to the Boulware Deser metric solution for Einstein Gauss-Bonnet theory with a 
cosmological constant [44-47]. Indeed setting a e // = 5/(3 as an effective Gauss-Bonnet coupling and A e ff = = 

—r)/f3 we obtain, 


/(r) = h(r) = k+ - r — *<Xeffk 2 eff + 16a r e / j • (3.12) 

where /z is the integration constant associated to mass. We see therefore that the effective higher dimensional coupling 
and the effective cosmological constant are rescaled via the definitions (3.9-3.10) by the scalar velocity parameter q. 
The above solution for the scalar field is valid for non planar horizons. For the particular case where k = 0, in order 
to have a solution we need to have that ij = and aplrj = 25<prj + f35A and then the scalar field has the form 

v>'(r) = ± ^q 2 + ww)- 



It is obvious that on the horizon of the black hole the radial component ip blows up but not the scalar field </>. 
Indeed in ingoing Eddington-Finkelstein coordinates 


the scalar field reads, 


u = t + 


f(r. 


■dr 


tp'(r) = ±- 


«±«v(i-«) 


(3.13) 


(3-14) 


We can now see that when f(r) vanishes the scalar field is regular in the future black hole event horizon (for the plus 
sign). A different sign choice renders the solution regular at the white hole horizon. 

If we now demand that the scalar field has no time dependence, in other words taking the limit q —> 0, we see that 
4>{r) = const. In order to examine the full spectrum of solutions for the scalar static case we must go back to the 
original equations and solve them again from scratch. 

We will again examine the case where f(r) = h(r). The only non vanishing component of J M is (3.3), and the 
scalar field equation admits a first integral. Following [36, 37] we can set the integration constant to zero and thus 
have a single equation which we solve for /(r), giving us again as a solution (3.8). We are now left with the metric 
field equations. We find two solutions, 


a = 

or 


a = 


2SCV + PSA ^ k = q andthus ISC-a/3 


Pv 


PS fir) 


c\ e 

iry^iSPC + 26A) and r\ = and thus <p'{ t) = 

op 


3/3 2 
~4S 


I SC-ap 

p5 f (r) 


(3.15) 

(3.16) 


As we see in both cases the solution for the scalar field is the same whereas the metric satisfies (3.8) for the relevant 
constants. Combining the constraints of the second case we see that we can derive the constraint of the first solution. 

o<o2 . . 

Still that case does not in any way imply the relation ry = -jj- and it is only valid for flat horizons, whereas the 

... . . qo2 ^ ^ ^ 

second solution is valid for any k. Still the relation ry = for non zero q is mandatory in order to have solutions 
for flat horizons. Additionally we see that when — ap = 0 the scalar field is constant and this is precisely the 
case mentioned already where q —> 0, since then the constraints (3.9),(3.10) produce the aforementioned relation. It 
is clear that when we abandon time-dependence, the scalar field blows up on the horizon of the black hole as was 
expected, [27]. 


IV. CONCLUSIONS 

In this paper we have found and studied spherically symmetric solutions of higher order scalar tensor theories. We 
have investigated not only 4 dimensional Horndeski theories but also theories including higher dimensional galileons 
stemming from derivative scalar couplings with the Lovelock tensor (1.3). We dubbed these Lovelock-galileon terms. 
With the help of these we have found a regular scalar tensor extension of the Boulware Deser solution [44-47]. The 
solution has the basic horizon structure of the afore mentioned solution but with a time and radially dependent scalar 
field which is furthermore regular at the black hole event horizon. As noted before, [27],this regularity feature is due 
to the linear time dependence of the scalar field. To investigate the role of time dependence of the scalar field further, 
we also studied a scalar field with only time dependence. In other words although our spacetime has a Killing time like 
vector the scalar field does not have this symmetry. We found spherically symmetric star solutions but they always 
have a singular event horizon, practically the scalar shifts the zero of one the metric functions but not the other thus 
producing a singular space time. Far from the event horizon radius solutions are GR like and asymptotically flat. 
Furthermore, as we let the scalar velocity charge go to zero we recover GR solutions with or without a cosmological 
constant. 

The scalar tensor version of the Boulware-Deser solution presents certain self tuning features reminiscent to those of 
de Sitter black hole solutions [27], [32]. Indeed we see that the bulk cosmological constant as well as the Gauss-Bonnet 
coupling are rescaled by the scalar velocity charge q. This could have interesting consequences in a Randall-Sundrum 
type braneworld scenario where the fine tuning of the brane tension with the bulk cosmological constant could be 
relaxed providing a 5 dimensional and regular self tuning solution. But in order to investigate such scenarios we would 
need to carefully study junction conditions and in particularly in the presence of higher order scalars terms. 
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Appendix A: Generalized deltas 


Here we give some definitions on generalized deltas. For more info look in [61] 


rCti CX.2 • ••CX. ri 

°/3i/3 2 .../3„ 




where in denoting the antisymmetry we include the normalizing factor, for example 


(Al) 


(A2) 


This way for example 


5 


OtlCX.2 

Pl&2 




(A3) 


Now making contractions we can define generalized deltas with fewer indices, but the dimensionality of the space that 
we are working on, plays a crucial role on the value of the numerical factors that appear, after the contraction has 
been made. So following [61] we have that 


raia2...asa a 4-i...a„ 

/3i/32-../3 s a s +i...a n 


(d-s)\ 

(d— ri)\ 


5 


a 1 a 2 ...a s 

Plfa—Pa ' 


Contracting fully the indices gives us 



ai -“" (d — n)\ ’ 


(A4) 


(A5) 


Appendix B: Lovelock Derivative Couplings 


It is easy to extend the derivative couplings for any Lovelock invariant. These type of couplings belong to the 
general type of higher order derivative couplings which can give second order field equations and extend Hordenski’s 
theory in arbitrary dimensions [49]. Each one of the Lovelock invariants can be written as 


C 


In ~ 



OL\OL2‘ • -OL2n-\Cx.2n 

Plp2"‘p2n-lp2n 


f?filp2 . . . f>(32n-lp2n 

L Oi 1&2 t O'2n-lOi2n 


(Bl) 


If we consider the above term in d = 2 n then this is a topological invariant and is not dynamical. For example for 
n = 2 we have (1.1), which is a topological invariant in four dimensions. Varying with respect to the metric gives us 
the following expression 


cl'Oil&2 “ ‘&2n— 1 C*2? 


TT V _ _ 

M ~ 2 n+l” nPl@2---P2n-102n ““1“2 


Z?/3i/3 2 

1 Irvt rvo 


, TD&2n — 1 /?2n 
L OL2n-lOi2n ' 


(B2) 


Then, if we consider the coupling, 


C 2 nH = 


(B3) 
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varying with respect to the metric gives us 

^ 2 ""„ = -\9v V HjV a <t>Vp<j> + H^V v) 4>V p <l> 

x\ct\CL2-- % Ol2n — lOt.2n R 02 "- 1 


2n+l 0 K/3lfi2--p2n-l(L-t' 

n 

2n+l ^K,0 1 0 2 ---02n-302n-2 02 


RPlPl ...^2n-3/3 2 n-2V^Vx(i 

il aiao il ao n _3a 9 v ■ v •*'r 


c Aai 0:2 • • - 0 : 211—3 <^ 2 n - 2 c* 2 n — 1 £* 2 ™ 

—«_ (ii9v)ol2 


v)ot2n-l&2n ''■^1^2 ^ —3 <2271 — 2 

P02nK Vpffi\(f) 

^V ft »V A 0. 


_rAaia2-<J2n-3«2n-2«2n-ia2n Dplp2 TJ02n-302n-2 y7K yy 

2n Kplp2---p2n-3p2n,-2 p2n ((1 S^Jain Jl aia 2 ' ' ' Jl a 2 „_ 3 a2„-2 V V 


(B4) 
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